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ABSTRACT 

We have recently shown ("see lBianco & Ruffin i 2004) that marked differences exist between the EQuiTempo- 
ral Surfaces (EQTSs) for the Gamma-Ray Burst (GRB) afterglows numerically computed by the full integration 
of the equations of motion and the ones found in the current literature expressed analytically on the grounds 
of various approximations. In this Letter the exact analytic expressions of the EQTSs are presented both in 
the case of fully radiative and adiabatic regimes. The new EQTS analytic solutions validate the numerical 
results obtained in Bianco & Ruffini (2004) and offer a powerful tool to analytically perform the estimates of 
the physical observables in GRB afterglows. 

Subject headings: gamma rays: bursts — Shock waves — Hydrodynamics — ISM: kinematics and dynamics 
— gamma rays: observations — relativity 



1. INTRODUCTION 
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We have recently shown (see iBianco & Ruffinil|2004 ll) that 
in Gamma-Ray Bursts (GRBs) marked differences exist be- 
tween the equitemporal surfaces (EQTSs) for the afterglow 
numerically computed by the full integration of the equa- 
tions of motion and the ones found in the current literature, 
expressed analytically on the grounds of various approxima- 
tions. Indeed, the approximate formulae in the current liter- 
ature overestimate the size of the EQTSs approximately by a 
factor of 27% or 20% in the early part of the afterglow (f^f = 35 
s), in the adiabatic and fully radiative cases respectively. Cor- 
respondingly, they underestimate the size of the EQTSs ap- 
proximately by a factor of 15% or 28% in the latest part of the 
afterglow (f^f = 4 days). The precise knowledge of the EQTSs 
is essential to obtain the observational properties of GRBs. 

In this Letter, progress is made in making manifest the dif- 
ference between the exact expressions of the EQTSs and the 
ones obtained by approximate methods in the current litera- 
ture: the exact analytic expressions of the EQTSs are found in 
the case of both fully radiative and adiabatic regimes. 
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dM = ^li^int + dMism 

^ dMism = A-Trmptiismr^dr 

where 7, and M are the pulse Lorentz gamma factor, in- 
ternal energy and mass-energy respectively, nism is the ISM 
number density, nip is the proton mass, e is the emitted frac- 
tion of the energy developed in the collision with the ISM and 
Mism is the amount of ISM mass swept up within the radius r. 
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where is the starting radius of the shock front. In general, 
an additional equation is needed in order to express the de- 
pendence of e on the radial coordinate. In the following, e is 
assumed to be constant and such an approximation appears to 
be correct in the GRB context. 



2. THE EQTS FOR GRBS 

2.1. The afterglow hydwdynamical equations 

The discovery of the afterglow (ICosta et al.ll 1997ft has of- 
fered a very powerful tool for the understanding of GRBs. 
Consensus has been reached that such an afterglow originates 
from a relativistic shock front propagating in the Interstel- 
lar Medium (ISM) and that its description can be obtained 
from energy and momentum conservation in relativistic hy- 
drodynamics. Consensus exists, as well, that the shock fluid 
is concentrated in a thin shell. The fulfillment of the en- 
ergy and momentum conservation in the laboratory reference 
frame leads to (see e.g. .Piraa.1999.: .Ruffini et aL.2003t and 
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2.2. The definition of the EQTSs 

For the case of a spherically symmetric expansion consid- 
ered in this Letter, the EQTSs are surfaces of revolution about 
the line of sight. The general expression for their profile, 
in the form d = i}{r), corresponding to an arrival time t„ of 
the photons at the detector, can be obtain ed from (see e.g. 
iRuffini et alJ2003tlBianco & Ruffinl200l : 

cta = ct(r)-rcos'd + r* , (3) 

where r* is the initial size of the expanding source, d is the 
angle between the radial expansion velocity of a point on its 
surface and the line of sight, and t = t(r) is its equation of mo- 
tion, expressed in the laboratory frame, obtained by the inte- 
gration of Eqs.Q. From the definition of the Lorentz gamma 
factor 7"^ = 1 - (dr/cdt)^, we have in fact: 

ct(r)= r[l-rHr')]-'''dr', (4) 
Jo 
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where 7(r) comes from the integration of Eqs.Q. 

It is appropriate to underHne a basic difference between the 
apparent superluminal velocity orthogonal to the line of sight, 
y-L ^ ^jjjj jjjg apparent superluminal velocity along the line 
of sight, v'l ~ 7^v. In the case of GRBs, this last one is the 
most relevant: for a Lorentz gamma factor 7 ~ 300 we have 
yll ~ lO^c. This is self-consistently veri fied in the structure o f 
the "prompt radiation" of GRBs, see e.g. iRuffiniet alJll2()0l . 

2.3. The case of an adiabatic regime 

We first examine the case of an adiabatic regime (e = 0). 
The dynamics of the system is described by the following 
solution of the Eqs.Q (see e.g. IPiraniri999t and references 
therein): 
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where 70 and Mg are respectively the values of the Lorentz 
gamma factor and of the mass of the accelerated baryons at 
the beginning of the afterglow phase. 

We haveperformed an exact analytic integration of Eq.(|^ 
using Eq.(|5) and, as a consequence, we have the exact analytic 
solution: 



t(r)= 7, 
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where to is the value of the time t at the beginning of the 
afterglow phase and m° = {A/3>)Timpni,^r\. 

The analytic expression for the EQTS in the adiabatic 
regime can then be obtained substituting t(r) from Eq.(|6j in 
Eq.(|3j- We obtain: 



COSl9 : 



(7) 



ctg ^ r 
r r 



(m°/MB 



2.4. The case of a fully radiative regime 

We turn now to the case of a fully radiative regime (e = 1). 
The dynamics of the system is given by the following solu tion 
of the Eqs.O (see e.g. iPiranlllQQQt iRuffini et alJl2003[ and 
references therein): 



7 = 



1 + (M.s,n/Mfl) (1+7;') [1 + (1/2) (Mi„n/Mfl)] 

To' + {Mi, J Mb) (1 + 7o') [1 + (1/2) {Mi, J Mb)] 



(8) 



Again, like in the adiabatic case, we have performed an ex- 
act analytic integration of Eq.(|4|i using Eq.ljSji. As a conse- 
quence, we have (Ruffini et al. 2003): 
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Fig. 1 . — Comparison between EQTS.S in the adiabatic regime (solid lines) 
and in the fully radiative regime (dashed lines). The upper plot shows the 
EQTSs for = 5 s, = 15 s, = 30 s and = 45 s, respectively from the 
inner to the outer one. The lower plot shows the EQTS at an arrival time of 2 
days. 



where A = ^{Mb- m°) /mf and C = Mb^Jo - l)/(7o + D- 

The analytic expression for the EQTS in the fully radiative 
regime can then be obtained substituting t{r) from Eq.(|9} in 
Eq.(|3j- We obtain: 
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2.5. Comparison between the two cases 

The two EQTSs are represented at selected values of 
the arrival time tg in Fig. [2 where the illustrative case of 
ORB 991216 has been used as a prototype. The initial con- 
ditions at the beginning of the afterglow era are in this case 
given by 70 = 310.131, ro = 1.943 x lO'"^ cm, fp = 6.481 x 10 ^ 
s, r* = 2.354 x 10** cm fsee lRuffini et alJ2001albLl2()()ll2()()l . 
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3. APPROXIMATIONS ADOPTED IN THE CURRENT LITERATURE 

In the current lite rature two different treatmen ts of the 
EQTSs exist: one by Panaitescu & Meszaros ( 1998) and one 
by Sari ( 199&) later appUed also by Granot et al. (J999) (see 
also Pira nll999tl2000Hvan Paradiis et alJ2000l and references 
therein). 

In both these treatments, instead of the more precise dy- 
namical equations given in Eqs.( l5l8> . the following simplified 
formula, based on the "ultrarelativistic" approximation, has 
been used: 

7cxr-", (11) 

where a = 3 in the fully radiative and a = 3/2 in the adiabatic 
cases. A critical analysis comparing and cont rasting our exact 
solutions with Ea. Jl 11 1 has been presented in Ruffini & Biancd 
(|2004). As a further approximation, instead of the exact 
Eq.0, they both use the following expansion at first order 
in 7"^: 



ct(r): 
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Correspondingly, instead of the exact Eq.(|6j and Eq.(|9}, they 
find: 
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The fi rst expression has been given bvl Panaitescu & MeszarosI 
([1.998) and applies both in the adiabatic (a = 3/2) and in the 
fully radiative (a = 3) case s (see their Eq.(2)). The second 
one has been given bv iSaril (^98) in the adiabatic case (see 
his Eq.(2)). Note that the first expression, in the case a = 3/2, 
does not coincide with the second one: ^ari ( 1998) uses a 
Lorentz gamma factor F of a shock front propagating in the 
expanding pulse, with F = \/27- Without entering into the 
relative merit of such differing approaches, we show that both 
of them lead to results very dif ferent from our exact solutions . 

I nstead of the exact Eas.(l3l. lPanaitescu & MeszarosI J1998I) 
and lSaril Jl998h both uses the following equation: 



eta = ct(r)- r COS d , 



(14) 



where the initial size r* has been neglected. The following 
approximate expressions for the EQTSs have been then pre- 
sented: 



1? = 2 arcsin 




27oCffl 



1 



cos ?? = 1 - 



167^ 



2a+l 

3" 



,(15a) 



(15b) 



The first expression has been given bv lPanaitescu & MeszarosI 
( 1998) and applies both in the adiabatic (a = 3/2) and in the 
fully radiative (a = 3) cases (see their Eq.(3)). The second 
expression, where = j{-d = 0) over the g iven EQTS and 
ri = 16^lcta, has been given bv ISaril (Il998h in the adiabatic 
case (see his Eq.(5)). 

In Bianco & Ruffini (2004) we have compared and con- 
trasted the results of the approximate expressions given in 
Eqs.( I15> with the ones based on the exact solutions, there nu- 
merically comput ed an d here given for the first time in ana- 
lytic form in Eas. (l7ll0> . 

4. CONCLUSIONS 

The formulae we have obtained are manifestly different 
from the ones in the current literature. They are valid for 
any value of the Lorentz gamma factor and they may be ap- 
plied, as well, to the physics and astrophysics of supernovae 
and active galactic nuclei. However, as suggested by the ref- 
eree, a word of caution is appropriate: the applicability of the 
thin shell approximation, used in deriving Eqs.(n), is likely 
to break down when the non relativistic Newtonian phase is 
approached. There, the swept up ISM mass is no longer con- 
centrated in a thin shell as exemplified, e.g., by th e Sedov- 
Taylor-Von Neumann solution (see e.g. Sedov 1993I)- 

The new EQTS analytic solutions validate the numerical 
results obtained in Bianco & Ruffini (2004). We have indeed 
verified the perfect agreement between the results of the nu- 
merical computations, presented there, and the new analytic 
results, presented here. 

From the numerical examples given in iBianco & Ruffinil 
(|2p04) it is also clear that differences exist between the cor- 
rect treatment and the approximate ones all along the GRB 
afterglow process: the approximate treatments systematically 
overestimate the size of the EQTSs in the early part of the 
afterglow and underestimate it in the latest part. 

The analytic results presented in this Letter, when ap- 
plied to a specific model of the shock front emission process 
( Ruffini et al. 2004) duly taking into account the ISM filamen- 
tary structure ( Ruffin i et al.»2005 ). allows to make precise pre- 
dictions of the luminosity in fixed energy bands and of the 
instantaneous as well as time integrated spectra of GRB after- 
glow. 



We are thankful to the anonymous referee for very good 
suggestions on the manuscript. 
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